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Abstract 



We provide a derivation from first principles of the primordial bispectrum of scalar per- 
turbations produced during inflation driven by a canonically normalized scalar field whose 
potential exhibits small sinusoidal modulations. A potential of this type has been derived in 
a class of string theory models of inflation based on axion monodromy. We use this model 
as a concrete example, but we present our derivations and results for a general slow-roll 
potential with superimposed modulations. We show analytically that a resonance between 
the oscillations of the background and the oscillations of the fluctuations is responsible for 
the production of an observably large non-Gaussian signal. We provide an explicit expres- 
sion for the shape of this resonant non-Gaussianity. We show that there is essentially no 
overlap between this shape and the local, equilateral, and orthogonal shapes, and we stress 
that resonant non-Gaussianity is not captured by the simplest version of the effective field 
theory of inflation. We hope our analytic expression will be useful to further observationally 
constrain this class of models. 
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The study of anisotropies in the cosmic microwave background radiation over the past two 
decades has dramatically improved our understanding of the early universe. There is now 
strong evidence that the anisotropies we see today originated from primordial fluctuations 
generated in the very early universe, and we have learned that these primordial fluctuations 
have a nearly scale-invariant spectrum. Furthermore, the data still contains no evidence 
for a deviation from adiabaticity, or Gaussianity pp. Even though the case is by no means 
closed, these properties of the primordial fluctuations certainly support the idea that they 
originated as quantum fluctuations during inflation, a phase of nearly exponential expansion 
of the universe [2]. 

While observations are now good enough to rule out some of the simplest inflationary 
models involving only a single slowly rolling field with canonical kinetic term, other models 
in this class are still compatible with all existing data. These models predict an adiabatic 
spectrum with primordial non-Gaussianities that are too small to be observed, but this is 
not a generic prediction of inflation. Many models of inflation have been constructed and 
studied that can lead to an observable departure from Gaussianity. Some popular possibilities 
are models with multiple fields, non-canonical kinetic terms, light spectator fields and a 
violation of slow-roll. Beyond an existence proof that observably large non-Gaussianities 
can be generated, these models provide us with useful theoretical expectations to guide our 
search. 
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For a Gaussian signal, all odd n-point functions vanish and the higher even n-point func- 
tions are given in terms of sums of products of the two-point function. The most straightfor- 
ward way to look for a departure from Gaussianity is then to look for a non-zero three-point 
function. In Fourier space the three-point function depends on three momenta. Translational 
invariance of the background geometry ensures that these momenta add up to zero and thus 
form a triangle. Rotational invariance furthermore dictates that the three-point function can 
only depend on the three independent scalar products of these momenta. The information 
contained in the three-point function can thus be captured by a function of three variables, 
that can be thought of as two angles and one side of the triangle. Since the dependence is a 
priori completely arbitrary, a model independent measurement would be desirable and would 
provide a precious criterion to discriminate between otherwise indistinguishable models. Un- 
fortunately, progress in this direction is very hard. (For a review see e.g. [3]). Essentially all 
phenomenological analyses start from some explicit form of the three-point function guided 
both by theoretical expectations and by the simplicity of the numerical analysis necessary 
to compare it with the data. (See, however, jl].) Once a "shape" has been chosen, only the 
amplitude of this type of non- Gaussianity remains as a parameter, which is conventionally 
called J sha P e . So far only a handful of scale-invariant shapes have been looked for in the 
data. The most recent observational bounds on the magnitude for various shapes from the 
7-year WMAP data at 95% CL are pQ: 

local non- Gaussianity — 10 < / local < 74, 

equilateral non-Gaussianity —214 < f C( i ml < 266 , 
orthogonal non-Gaussianity —410 < / ortho < 6 . 

Planck data will make it possible to tighten the error bars by about a factor of five, and we 
may soon find out whether it is necessary to go beyond the simplest models of inflation. 

As already briefly mentioned, departures from the slow-roll condition can potentially 
lead to large non-Gaussianities. Two conceptually distinct possibilities were first explored 
by Chen, Easther, and Lim. An otherwise smooth potential might either exhibit a sharp, 
localized feature [5] (see also [6]), or it might display a periodic modulation that averages 
to zero [7]. Chen, Easther and Lim have performed a numerical analysis of both scenar- 
ios 07]. They show that a large non-Gaussian signal can be produced without violating the 
constraints on these models from measurements of the two-point function. They also pro- 
vide a heuristic estimate of the signal for equilateral configurations for the case of resonant 
production. 

In this work, we analytically compute the scalar primordial bispectrum generated for a 
modulated potential for arbitrary momentum configurations from first principles. Our work 
is motivated by a class of models derived from string theory [8119] . but let us stress that 
these are not the only models in which such oscillations are expected to arise. In large 
field inflation, the inflaton potential must be flat over a range in field space large compared 
to Mpi. From the point of view of effective field theory, this seems unnatural unless there 
is an underlying shift symmetry. Axions are thus natural candidates for the inflaton in 
large field inflation. It then seems plausible that the potential might receive small periodic 
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contributions from non-perturbative effects. These periodic contributions might be due to 
instantons in a gauge sector the axion couples to, or, in the context of string theory, they 
might arise from Euclidean branes or world-sheet instantons. Whether string inspired or not, 
as soon as we invoke the shift symmetry of axions to explain why the inflaton potential is 
so flat, we should admit the possibility of small periodic modulations in the potential which 
may lead to observational consequences. To be specific, the potentials we will consider are 
of the form 



V{<j>) = V o (0)+A 4 cos 



./ 



where is a canonically normalized real scalar field, and Vo(<p) is assumed to admit slow-roll 
inflation in the absence of modulations, i.e. for A = 0. A particular model of this form with 
Vo(0) = /i 3 was obtained from a string theory construction in [8ll9] , and we will sometimes 
focus on this special case for concreteness. The parameters A and / have dimensions of a 
mass and are a priori undetermined. However, consistency of a more fundamental description 
of the system, in our case string theory, will typically limit them to lie in a certain range. We 
will assume that the potential is monotonic at least near the values of the scalar field around 
which the modes we observe in the cosmic microwave background (CMB) exit the horizon 
and do not consider models in which the inflaton gets trapped. This can be summarized by 
requiring that the monotonicity parameter 

A 4 

Except for a linear potential, this parameter depends on the value of the field. We will 
evaluate it at <p = the value of the scalar field at the time when the mode with comoving 
momentum equal to the pivot scale k = k* exits the horizon. To be compatible with WMAP 
data, the monotonicity parameter must satisfy 6* 1 for both linear [H] and quadratic 
Vq{4>) [10] • Other potentials of this form have not been compared to the data, but we expect 
this to be true for the general case and treat 6* as an expansion parameter. 

Our main result is that the three-point function of scalar curvature perturbations to 
linear order in 6* at some late time t, when the modes have exited the horizon, takes the 
form0 

1 



(7e(k 1 ,t)^(k 2 ,t)^(k 3 ,t)} = (2tt)'A 



sm 



f 



with 



k-^ k 
hiK/K ) + 



2 p 
2^3 



^(kx+ka + ks) 



Eh 



T 



In Kj K + 



36, V2n 



/ 2e. 



3/2 



;i.3) 



;i.4) 



1 We have set Mpi = 1 in this formula and will do so throughout the paper. As usual, it can be re-inserted 
by dimensional analysis. 
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where e* denotes the value of the slow-roll parameter derived from the smooth part of the 
potential Vq(4>), evaluated at the time the pivot scale fc* exits the horizon. The quantities 
and e* are model dependent, but are easy to calculate for any given model. The comoving 
momentum K = ki + k^ + k^ is the perimeter of the triangle in momentum space. The dots 
in (II. 31) stand for terms that can be neglected either because they are suppressed by higher 
powers in the slow-roll parameters for the smooth part of the potential or by positive powers 
of //i/2e*. From (11.41) one can see that large non-Gaussianity requires / /y/2e^ «C 1. So this 
will be the regime of interest in this the paper. The second term is suppressed compared to 
the first by a factor of f / y/2e 1f . It is negligible except for squeezed triangles where one of the 
momenta is much less than the other two. It ensures that the consistency relation of [TT] 
(see also [12]) holds. 

We compare our analytic result with the numerical analysis of [T3] for the linear potential 
of axion monodromy. We find agreement with their numerical results for the bispectrum at 
the per cent level. We compare our results with the numerical analysis of [7] for the quadratic 
potential. Identifying their parameter P 2 with 9A^(/c*)/10, we find agreement with their 
numerical results as well. 

This type of non-Gaussianity, which we refer to as resonant non-Gaussianity following 
the nomenclature in [7], is nearly orthogonal to all commonly studied shapes. The cosine 
defined in [14115] is less than 10% for the entire range of parameters relevant for axion 
monodromy inflation. This is intuitively clear. The shape of resonant non-Gaussianity 
rapidly oscillates around zero while the other shapes are slowly varying. As the number 
of oscillations increases, the cosine decreases. The number of oscillations over the scales 
observed in the cosmic microwave background is approximately given by \/2e*/ /. For our 
string theory example, this implies that an axion decay constant / = 10 -3 leads to about 
90 periods in the cosmic microwave background. The cosine with local, equilateral, and 
orthogonal shapes is then around 2%, and the largest amount of non-Gaussianity for this 
value of / that is consistent with the observational constraints on the two-point function is 
/ res ss 80. 

As we will explain, the resonant effect we are studying arises from a term in the interaction 
Hamiltonian that is higher order in the slow-roll parameters. As a consequence, it is not 
captured by the simplest version of the effective field theory of inflation [16] (see also [T7]). 
A more detailed discussion is postponed to Section [5j 

The paper is organized as follows. In Section [21 we derive the time evolution of the 
curvature perturbation for a potential (11.11) and obtain the power spectrum. The derivation 
is independent of the ones presented in [9]. It agrees with the results that were found there 
for the linear potential. In Section [3j we calculate the bispectrum and check that it fulfills 
the consistency relation of [TT] for squeezed triangles. In Section HJ we compute the overlap 
between the shape of resonant non-Gaussianity and the three most common shapes in the 
literature that have been compared with the data. We show that it is very small for the 
range of parameters relevant for our string theory example. Section [5] contains a discussion 
of our results. In Appendix |A], we present the details of the calculation of the inflationary 
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background solutions for the potential (11. ip . Appendix [B] provides the relation between 
various different definitions for slow-roll parameters commonly used in the literature. 



2. The Mode Functions and the Power Spectrum 

As we show in Appendix |Aj the background solution for a scalar field with potential (II. ip 
can be derived to first order in 6* and to leading order in the slow-roll parameters of Vq. In 
the limit / <C \f2^, which, as we shall see, is the regime where large non-Gaussianities are 
generated, the solution is well approximated by 

(j){t) = cj) (t) = sin — — . (2.1) 



^ \ f 

Here 0* is the value of the scalar field when the pivot scale k = k* exits the horizon^ e* is 
the slow-roll parameter in the absence of modulations evaluated at <pQ — <P*i an d 0o is the 
solution for the scalar field in the absence of modulations, i.e. for 6* = 0. It can be obtained 
as a function of time by integrating its equation of motion, but we will not need the result 
at this time. 

Since we will need them later, let us give the expressions for the Hubble slow-roll param- 
eters§ 

e = -^ and 5 = -^-. (2.2) 
H 2 2HH 

For the potential (ll.ip . it is convenient to calculate them in an expansion in the parameter 

b,, 

e = eo + ei + (bl) , (2.3) 
6 = 6 + 6 1 + O{%). (2.4) 

In the slow- roll approximation for O , and for / <^ \/2e*, one finds 

eo = e*, 5o = e* - 77* , (2.5) 
ex = -36,/v^cos , (2.6) 



/ 

^ = -36,sin(^ , (2.7) 

where 77* and e* are the values of the potential slow-roll parameters r] Vo = Vq /Vq and ey = 
(Vq/Vq) 2 /2 derived from the smooth part of the potential Vq, evaluated at the time the pivot 



2 For numerical calculations, we will take fc* = 0.002 Mpc -1 . The value of cj>* is model dependent. For the 
linear potential we use </>» = 10.88. 

3 For the convenience of the reader we have collected in other possible definitions of the slow- roll parameters 
and formulae for the conversion in Appendix [B] 
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scale k* exits the horizon. Notice that in this regime both e\ <C 1 and 8\ <C 1 as long as 
6* <C 1. On the other hand this is not the case for higher slow-roll parameters. For instance, 
one has 

8 X /H = 3b^ca*(^\ , (2.8) 

which becomes large for small / / y/2e^. 

Let us now turn to the spectrum of scalar perturbations. As in [9], we will choose a slicing 
such that 5<fi(x., t) = 0, and use a spatial diffeomorphism to bring the scalar perturbations in 
the spatial part of the metric into the form 

6 gij (x,t) = 2a(t) 2 H(x,t)6 ij . (2.9) 

The translational invariance of the background makes it convenient to look for the solution 
superposition of Fourier modes 

—TlQtrfe*"* , (2.10) 

where x are the comoving coordinates and k denotes the comoving momentum of the mode. 
Rotational invariance together with reality (or Hermiticity) of TZ(x, t) imply that the Fourier 
components of the most general solution can be written in the form 

K(k,t) =n k (t)a(k)+ni(t)a\-k), (2.11) 

where k is the magnitude of the comoving momentum k, a(k) can be thought of as a 
stochastic parameter in the classical theory or as an annihilation operator in the quantum 
theory, and TZk{t) is the mode function. When thought of as creation and annihilation 
operators a^(k') and a(k) satisfy the commutation relation 

[o(k), a t (k')] = (2vr) 3 5 3 (k - k') . (2.12) 

The time evolution of the mode function TZk{t) is governed by the Mukhanov- Sasaki equa- 
tion 118,191. For small e, it can be written in the form |20] 



gfr_2(l + fc + *)jK* +R> = . 
ax 2 x ax 

The initial conditions are such that for x ^> 1 

K k (x) ^ -JL^e* , (2.14) 
V 2ka(p 

where we have used the notation x = ~kr with conformal time r defined as r = f -M^. 
The Mukhanov- Sasaki equation implies that for i C 1 the mode function TZk{x) approaches 
a constant which we denote by TZ^\ where the superscript (o) indicates that the mode 
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is outside the horizon. It is related to the quantity A^(/c) that is commonly quoted to 
parameterize the primordial scalar power spectrum by 



(o) 



In the slow-roll approximation, i.e. for e < 1, (5 < 1, and assuming 5/H is small compared 
to both e and 5 

H 2 (t ) 

AiU&) = — 7—. — r (slow-roll approximation) , (2-16) 
8ir 2 e(t k ) 

where t k is the time at which the mode with comoving momentum k exits the horizon. 
However, as already pointed out in [9], the slow- roll approximation breaks down in models 
with modulated potentials because the magnitude of 5/H is no longer of quadratic order 
in e and 5. Furthermore, the slow-roll parameters are oscillatory functions whose frequency 
changes in time. When the frequency of this oscillation passes through twice the natural 
frequency of the mode, which is set by the momentum of the mode, parametric resonance 
occurs, which is not captured in the slow-roll approximation. 

In [9] our main concern was the power spectrum so that only the asymptotic behavior of 
the mode function was needed but not its detailed behavior as a function of time. The main 
concern of this work is the calculation of the bispectrum for which the knowledge of the 
time dependence of the mode functions is important. So let us calculate it. We will neglect 
the effect of e and 5 in equation (12 . 1 3[) for simplicity, but they could be restored without 
too much extra trouble. Furthermore, we will make use of the fact that the amplitude of 
e\ is suppressed compared to that of 5\ by a factor of fy/2e^ and drop it as well. The 
Mukhanov-Sasaki equation ( 12 . 1 3[) then becomes 

d 2 n k 2(1 + 5 1 (x)) dn k 

—3—5 1 r Kk — U . [Z.L( ) 

dx A x dx 

As was shown in [9], for the linear potential parametric resonance occurs around x res = 
1/(2/0*). As we will see, for the general potential it happens at x res = y / 2e^/(2/). For 
x x res , i.e. much before the resonance occurred, we know that the effect of 5\ is negligible. 
Therefore the solution i^| 

K k {x)=K^i^/*HV{x), (2.18) 

where 1Z k °\ is the value of 7Z k (x) outside the horizon in the absence of modulations and is 
fixed by the initial condition ( I2.14I) . and 



i x l^x 3 / 2 H^(x) = (l-ix)e ix . (2.19) 



4 To keep the dependence of the mode function on the slow- roll parameters eo and 5q, one should replace 
3/2 by 3/2 + 2cq + 6q in the discussion below. 
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Similarly, for x <C x res , i.e. long after the resonance has occurred, the background frequency 
is too high for the mode to keep up with it, and the effect of 8\ is again negligible. The 
solution there must take the form 



TZk(x) 



(°) 

fe,0 



,(+), 



'3/2 



'3/2* 



A slight generalization of our derivation in [9] implies that at late times cjj. = 1 
then seems natural to look for a solution of the form 



K k (x) = n 



(o) 
fe,0 



Ul -x^H^Jx) - <£\xM J-x^Hfhx) 



3/2* 



3/2* 



(2.20) 
0{bl). It 

(2.21) 



where c\r\x) vanishes at early times and goes to c£~ at late times. The Mukhanov- Sasaki 
equation then turns into an equation governing the time evolution of cj~\x). To linear order 
in this equation is 



d 

dx 



-2ix 



d 



1 — 7T '■■ ' { ' r] 



+ e 



2ix * d (-)/ \ 

x*dx k [ } 



-2i 



Si(x) 



X 



(2.22) 



For large x, which is where the resonance occurs as long as // y/2e^ <C 1, equation (I2.22p can 
be written in the form 



d_ 

dx 



6 dx Ck [X) 



-2i 



8i(x) 



X 



Using the expression for 8\{x) given by equation (12. 7ft together witrjfl 

4>o(x) = 4>k + v / 2eTlna; where 4> k = 0* - v / 2e7hi k/k* , 



(2.23) 



(2.24) 



where is the value of the scalar field when the mode with comoving momentum k exits 
the horizon, this can immediately be integrated once to give 



-6ib-)=e 2ix 



< os ( — H — inx 



(2.25) 



We will be able to use this result in the next section to argue that the modification of the 
mode function can be ignored when calculating the bispectrum to leading order in 6*. Let 
us now integrate this equation once again to get a better idea for what the function c k \x) 
looks like. To separate the leading and subleading contributions, it is convenient to write 

c ( ~\x) as 



c { -\x) 



-3ib- 



f 



e t 



i I 2ix'-i^\nx' , A 

dx e f + e f 



i I 2ix'+i^±-\nx' 

dx e i 



(2.26) 



After some manipulations, these integrals can be recognized as incomplete T-functions or 



5 See Appendix [Al for a details. 
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Figure 1: The top left shows the dominant contribution to Re[c^, (x)] coming from the 
first integral in equation (12. 26ft . while the top right shows the subdominant contribution 
to Re[cjj. \x)}, which comes from the second integral in equation (12. 26ft . The bottom left 
shows their superposition. The bottom right shows the superposition as black dotted line 
compared to the numerical solution in orange. As can be seen, the two results are essentially 
indistinguishable. The dashed line represents the asymptotic value for x going to zero (i.e. 
outside the horizon) found by performing the integral in f!2.28j) in the stationary phase 
approximation. All plots are for a linear potential with / = 10~ 3 , b = 1CT 2 , 0* = 10.88, and 
for a value of comoving momentum such that 4>k = 10.7. Notice that Re[c£ \x)] changes 
around x ies = 1/2/0* « 45 as expected. 



closely related exponential integrals. The result can be written as 



,(-), 



x 



2 x/27: 



e f 1 



T+i In 2 



1 f r 



/ 2e^ 



f ' 



1 + 



-2ix 



f 



+ 



—2ix 



(2.27) 



To gain some intuition, the stationary phase approximation is useful. The phase of the first 
integrand in (I2.26P will become stationary near the resonance at x — x res , while the phase of 
the second integrand is never stationary (since x > 0). The first term on the right hand side 
of equations (12.26jl and (12.27P will thus be the dominant contribution, with the subdominant 
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contribution from the second term exponentially decreasing for decreasing //\/2e*. Figured] 
shows the leading and subleading contributions to the real part of c£~ (x), their superposition, 
as well as the comparison to the numerical result for the c* k \x) from its evolution equation 
( 12.22)) for a linear potential. Equation f |2.22j) is valid for arbitrary x and thus exact at linear 



order in 6*. As a byproduct, we have given an expression for c 



(-) 

k 



,(-) 



(0) of the form 



•, :>ub-L< i: fdxe** ' ' 



(2.28) 



When evaluated exactly as in (I2.27P or using the stationary phase approximation, one finds 
that up to an unimportant phase 



1/2 



e f . 



(2.29) 



For the linear potential we can set ^f2el_= 1/0*, and we find that this expression agrees 
with equation (3.40) in [9 J for small f4>*u Using equations (I2.2ip and f |2.29|) as well as the 
behavior for the Hankel functions for small arguments, one finds that the primordial power 
spectrum of scalar fluctuations is of the form 



\n 



(°)|2 



\n 



(°)|2 
fc,0l 



1 + 36* 



2nf_ 



1/2 



COS 



/ 



(2.30) 



or equivalently 



A 2 



K 



(k) 



A 2 



n 



{k* 



1 + 5n s cos 



./ 



with 5n s = 36* 



2tt/ 



1/2 



(2.31) 



where once again is the value of the scalar field at which the mode with comoving mo- 
mentum k exits the horizon, e* is the value of ey when the pivot scale k = k* exits the 
horizon, and 6* = A 4 /V^(0*)/. We have restored the dependence on e and 8 through the 
appearance of the scalar spectral index n s , which in the approximation we are using is given 
by n s = 1 — 4e — 25 = 1 — 6e* + 2r/*. 

Everything we have said here about the primordial power spectrum for the scalar modes 
is valid for small //\/2e*, which, as we will see, is the regime in which observable non- 
Gaussianity can be generated. In [9], the interested reader can find the result for the linear 
potential for general /</>*. 

We have performed numerical calculations to check these analytic results, and we find 
good agreement. Most of the numerical calculations were done for the linear potential 
relevant for axion monodromy inflation, but we have also performed some checks for the 
case of a quadratic potential as well as Vo(0) = /i lo / 3 2//3 motivated by [2T] . Our numerical 
results for the amplitude of the modulations as well as the frequency agree with our analytic 



3 This is true after dropping the same independent phase that was dropped there. 
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result at the per cent level in all cases. The discrepancy between our analytic result and the 
numerical results for the power spectrum in [13] can be traced to an initial value for k/aH 
in their numerical calculation that was too small to capture the resonance for small axion 
decay constants. This issue will be easy to fix. 



3. The Bispectrum 

Let us now turn to the calculation of the three-point function. To leading order in pertur- 
bation theory, the three-point function in the "in-in" formalism [22] (see also [23p4] and 
references therein) is given by 

(7e(ki,t)^(k a ,^(ka,t)) = -i / ^([^(ki^^ka^^ka,*),^)]), (3.1) 

J — oo 

where the expectation value is taken in the in-vacuum, and the interaction Hamiltonian 
Hi was first worked out in [11] (see also [25|26] ). The term responsible for the dominant 
contribution in models with oscillatory potentials in our notation is given by [7] 

flj(t) D - J d 3 x a 3 (t)e(t)5(t)ft 2 (x, t)TZ(x, t) . (3.2) 

Using equations (12.101) and ( 12. lip , one finds that the contribution to the three-point function 
from the term (13.21) in the interaction Hamiltonian is given by 

(ft(k!, t)K(k 2 , t)K(k 3l t)) = (2tt)H S 3 (k t + k 2 + k 3 )U kl (t)K k2 {t)K k , (t) x 

I dt' 2a 3 {t')e{t')S{t') \n* ki {t')K* k2 {t')TZ* k3 (f) + 2 perm] + c.c. (3.3) 

J — oo 

Observational constraints on the two-point function imply that knowing the result to linear 
order in 5* will be enough. The oscillatory nature of S is what makes this contribution 
the dominant one. So to linear order in &*, we can replace e by eo, 5 by Si and use the 
unperturbed mode functions for lZ k (t), i.e. equation (12.181) . One might be concerned that 
the derivative of the correction to the mode function becomes large during the resonance and 
should be kept, but this is not the case. To see this, notice that there is no contribution to 
the integral after the modes have frozen out. In fact, we will see that the main contribution 
arises when the modes are still deep inside the horizon. In this limit, i.e. for large x, the 
ratio of the absolute value of the time derivative of the unperturbed part lZ kj o(t) and the 
absolute value of the time derivative of the correction 1Z kt i(t) becomes 



\n 



k,l\ 



(—)/ \ ■ d (-) / \ 
k ^ +% ~dx~ k ^ 



(3.4) 



The results in the last section imply that this is small. Equation ( 12. 27ft tells us that the abso- 
lute value of the first term is never significantly larger than 36* ^/7r/2(/ /a/26*) 1 / 2 , and equa- 
tion ( I2.25P reveals that the absolute value of the second term is always less than 6bf /y/2e*. 
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As we will see, large non-Gaussianities can only be generated for decay constants satisfying 
/« 



x 2e*, so that both are small. We conclude that 



\n 



k,l\ 



\K 



«1. 



(3.5) 



fe,0 



In other words, we can use the unperturbed mode functions, and the approximation becomes 
better for decreasing axion decay constant. We are interested in the value of the three-point 
function after horizon exit. In this case, we can replace the factors TZ^if) outside the integral 
by TVj^o, and take the upper limit of the integral to zero. We will drop the dependence of the 
mode functions on the slow-roll parameters e and 8 , as well as e\, and use expression (I2.18P 
for the mode functions inside the integral. This leads to 



(ft(kx,*)ft(k a ,*)rc(k3,*)) = (2vr) 7 A 



7 A 4 1 

ft ^3l3l3 



' U-$ U'J h"^ 
rv-^ rL<2 rue 



<5 3 (ki + k 2 + k 3 



(3.6) 



dX-^-e' iX 
8H 



-ik x k 2 k 3 -jYl k i k i + -j^ K ( k i + fc 2 + fef) 



+ c.c . 



where H and 5\ should be thought of as functions of X = —Kt, and K = k\ + k 2 + k 3 is the 
perimeter of the triangle in momentum space. 

When visualizing the results, it is often convenient to introduce a quantity that contains 
the information about the deviation of the three-point function from scale invariance rather 
than the three-point function itself. We will use the notation of [7] and defind^ 



(Jl(k 1 ,t)K(k 2 ,t)'ll(k 3 ,t)) 



g(fcl,fc 2 ,fc 3 ) 

hk 2 k 3 



(3.7) 



It can then be seen from equation (13.61) that for all models whose three-point function receives 
its dominant contribution from the term in the interaction Hamiltonian f 1 3 . 2 j) . and for which 
the mode functions are well approximated by the unperturbed ones, one has 



Gjh, k 2 , k 3 
k\k 2 k 3 



i 

dX^e~ iX 
H 



-i — 



Y ^ 



X ^ ~kj + 'X 2 



K(k\ + k 2 2 + k 2 3 ) 
k\k 2 k 3 



+ c.c. (3.8) 



If the integral receives its main contribution from a small neighborhood around some 
value of X = X res , as is the case in our example, we can replace 1/X and 1/X 2 by l/X Tes 
and 1/X 2 es , respectively, and find that the shape of non-Gaussianity is given by 



G(h, k 2 , k 3 
k\k 2 k 3 



ImX, 



K 



1 v - ^ ki 1 K{kj + k 2 + k 2 

-2^TT ReI K-Y-2 h k 2 k 3 



-ImX 



A 



(3.9) 



7 There is a factor of 9/10 between our definition of Q and theirs, i.e. <7therc = 105hcrc/9, which was 
introduced there presumably to match the WMAP conventions for the local case, where a famous factor of 
3/10 appears. The remaining factor of 3 arises because the matching is conventionally done in the equilateral 
limit where three terms become equal. 
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where we have defined the integral 

Ik 



dX%e~ iX 
H 



(3.10) 



We see that all we need in order to calculate the shape of non-Gaussianities for models 
of this class is the quantity 5\/H as a function of X. We have already given this quantity 
as a function of the scalar field in equation (12. 8p . It remains to write it as a function of X. 
To leading order in the slow-roll approximation the scalar field is given in terms of X as 

4> (X) = (j) K + y/2e^\nX with (j) K = 0* — y/2el In K/k* , (3-H) 



where <pK is the value of the scalar field at the time the mode with comoving momentum K 
exits the horizonjf) This gives 



Si 
H 



cos 



+ 



InX 



(3.12) 



/ \f f 

The three terms in the integral (13. 8 j) can then be recognized as T-functions, and the inte- 
gral (13. 8p can be done analytically. In the regime / <C \/2e*, in which the resonance occurs 
deep inside the horizon, we can also use ( 13.91) . The integral (13. lOf) then takes the form 

oo 



1, 



K 



f 



dX e 



cos 



which can be written in terms of T-functions 



1, 



K 



V 



e~^f~Y I 1 + i 



/ 



e / 



9K 

f 



InX 



+ e 1 



./ 



iiK. 



(3.13) 



(3.14) 



The absolute values of the T-functions in the first and second term are identical so that the 
first term dominates for small // \/2e* because of the exponential factors. This dominant 
contribution arises from a neighborhood of size (y/2e^/ f) 1 ^ 2 around X res = y/2e*/ f where 
the phase of the integrand in equation (13 . 13[) becomes stationary. Equation (13. 9 p can then 
be used as long as / <C y/2e*, which is the regime we are interested in. Either performing 
the integral directly in the stationary phase approximation or using Stirling's approximation 
in equation (I3.14p . one finds that up to a K- independent phase 

3/2 

1k = -~-^-[^\ e'-f. (3.15) 



2 V / 

Combining this with equation (13. 9p . we see that the shape of resonant non-Gaussianity is 
given by 



Gjh, fa, h 
hk 2 k 3 



3V27T&* 



/ 



.3/2 



sin 



4>K 

f 

f 
x/2Z 



f 



f2l 



■a. ■ 3 



Vk_ 
f 



K(k\ + k\ 



sin 



f 



(3.16) 



3 See Appendix [Al for details. 
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Other terms in the interaction Hamiltonian also contribute at order (/ /y/2e*) 2 , but these 
contributions are too small to be phenomenologically interesting and we will drop them. 
Ignoring a /^-independent phase, we thus write the final result for the resonant shape as 



G(fa, fa, fa) 
fafafa 



f 



sm 



/2e» 
/ 



In Kj fa 



f 



or equivalently for the three-point function 



(3.17) 



(ftCki, t)ft(k 2 , t)n(ks, t)) = (2n) r A 4 nm ^ 5 3 (k! + k 2 + k 3 ) 

fa fc 2 % 



sm 



I2e* 
f 



In Kj fa 



with 



36* ( JTe 



f 



3/2 



(3.18) 



(3.19) 



The dots stand for terms that have been dropped because they are higher order in slow-roll 
or f/V2T*- 

A few comments on this result are in order. Notice that both the frequency of the 
oscillation and, when written in terms of the monotonicity parameter the amplitude / res 
depend only on y/2e*J f. If this type of non-Gaussianity were measured, it would thus not be 
possible to distinguish between different potentials from this measurement alone. However, 
a measurement of the amplitude of tensor modes would give us a direct measurement of e* 
and hence break the degeneracy!^] 

Concerning detectability, notice also that when the axion decay constant becomes too 
small, the frequency of the oscillations becomes too high to be experimentally resolvable. 
Let us assume the signal can be resolved if the period is longer than A£ ~ 1 around the 
first Doppler peak, i.e. near £ ~ 2000 The periodicity near a given value of £ is given by 
2n£f I \/2e*. So our condition is y/2e^/f < 2n£, and taking £ = 200 leads to 



f es <4x 10 4 6*. 



(3.20) 



In order to go from this to a numerical value for the upper bound on / res , one needs an 
upper bound on 6* from comparison of the predicted power spectrum with the data. So far 
this comparison has only been done for the linear potential [9], which will be the subject 



9 For a general discussion of the implications of a measurement of tensor modes for the inflationary theory 
see [27] and references therein. 

10 This will then also be true for all larger I 
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of the next subsection, and the quadratic potential [10] However, in [10], only decay 
constants / (/3Mpi in their notation) larger than 5 x 10 -3 were considered, which is too large 
to give a sizable / res . 

It may also be interesting to consider a potential of the form 



V(</>) = V (<f>) 



1 + A cos 



./ 



(3.21) 



In the approximation we have been working in, our results can immediately be translated to 
this type of potential by replacing b*fy/2e^ — > A. One finds that@ 



3AV27T 



1/2 



(3.22) 



This type of potential with V ((p) = |m 2 2 was studied by Chen, Easther, and Lim in [7]. 
One has 



/ 2^= — , 



(3.23) 



so that our result is 



3v^F A 



4 PVW* 



(3.24) 



This is smaller than their analytic estimate for the amplitude in the equilateral limit by a 
factor of 10a/7t/27, which is in agreement with their statement that their analytic estimate 
overpredicts the numerical result by about 30% @ 



3.1. The bispectrum for a linear potential 

We will now consider axion monodromy inflation [81:9] in some detail as a special case. The 
low energy effective theory describing the system is that of a canonically normalized real 
scalar field with potential 



V(<f)) = /i 3 + A 4 cos 



/ 



+ bf cos 



./ 



(3.25) 



The parameter /i has dimensions of a mass and is fixed by COBE normalization to be 
approximately /i ~ 6 x 10 -4 . 



11 For another observational constraint on an oscillatory power spectrum due to a deviation from the 
Bunch-Davies vacuum see [28] . 
12 See Appendix [Al for details. 

13 For the comparison, recall that there is difference by a factor of 10/9 between their definition of / rcs and 

Our / rcs from £ t herc = 10£horo/9. 
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For this potential, one finds = 1/0*, and 6* = b is now independent of 0*. We thus 
conclude that the shape of non-Gaussianities in this model is 



g(fci,fc 2 ,fc 3 ) 
hk 2 k 3 



f res sin 




or equivalently for the three-point function 



(^(k 1 ,t)^(k 2) t)^(k 3 ,*)} = (2vr) 7 A 4 



h 2 h 2 h 2 



5 3 (k! + k 2 + k 3 ) 



x/ 



■res 



sm 




where 



■res 



(3.28) 



8(/0*) 3 / 2 ' 



Once again, this result is valid at leading order in an expansion in 6, which is constrained to 
be much less than one by comparison of the two point function with CMB data [9] (pK is 
the value of the field at which modes with comoving momentum K = k± + k 2 + k% exit the 
horizon, 0* is the value of the inflaton field when the pivot scale fc* exits the horizon. Provided 
this happens around 60 e- folds before the end of inflation, one has 0* ~ 11. Finally the axion 
decay constant / is a free parameter. A typical value in an explicit string construction is 
10 ~ 4 < / < 10 _1 . The result is valid provided /0* <C 1. Logarithmic dependences on k 
which arise from the dependence of the mode functions on eo as well as 5q were neglected. 

As the axion decay constant decreases, the frequency of the oscillations increases linearly 
and, keeping the parameter A in the potential f!3.25j) fixed, the amplitude increases rapidly 
like / 5 / 2 . It is a natural question to ask for what values of the axion decay constant an 
observably large signal can be generated in this model while satisfying the bounds on the 
power spectrum from the data. As long as / < 10~ 2 , the bound at 95% confidence level is 
summarized approximately by bf < 10~ 4 [9\. Combining this with equation (13.281) . one finds 
that 



The regime in which the model can simultaneously be consistent with the constraints on the 
two-point function and generate an observably large three-point function is thus / < 6 x 10~ 3 . 
As we mentioned for the general case, requiring that the period of the oscillation should be 
larger than A£ ~ 1 for £ pa 200 leads to a lower bound on /. For the linear case, it is 
/ > 10~ 4 . The range of axion decay constants for which observably large non-Gaussianities 
can be generated is thus approximately 




(3.29) 



10 



4 </<6x 10~ 3 . 



(3.30) 



This is true provided the potential is monotonic, i.e. b < 1. 
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Notice that in this range <C 1 is always satisfied. 

The shape of resonant non-Gaussianity for axion monodromy inflation is shown in Fig- 
ure[2]for b = KT 2 , /0* = 2 x 10~ 2 , and fixed k x = K = 0.002 Mpc -1 . We chose this value of 
/ because both the leading contribution and the subleading contribution in are clearly 
visible. Notice that as the value of k\ changes, the phase of the oscillation changes. 



f0 t =O.O2M P 2 1 




Figure 2: This plot shows the shape G(ki, k 2 , k^)/{kik 2 kz) of resonant non-Gaussianity for 
the linear potential of axion monodromy inflation with b = 10~ 2 , = 2 x 10~ 2 and fixed 
k± = k* = 0.002 Mpc -1 . We use the notation x 2 = k 2 /k\ and X3 = k^/k\. The triangle 
inequality implies x 2 + x% < 1 and the quantity is symmetric under interchange of x 2 and X3 
so that we show in the plot only the region 1/2 < x 2 < 1 . 



We find that our analytic result forf res agrees with the values obtained by numerical 
integration in [13] at the per cent level] 15 ! 

3.2. Consistency relation 

As pointed out in [11] (see also p2]), in the limit in which one of the momenta, say, k% 
is much less than the other two, which are then roughly equal, fc 3 C ^ « fc 2 = k, the 
three-point function is related to the two-point function by a consistency relation 

lim ^ki.t)^,*)^,*)) ~ -i^l 2 -^— -^(k^T^W)}, (3.31) 
^3^0 J ti [tkj atk 
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For the comparison, notice that j!3| uses a momentum dependent quantity Jnl- In the equilateral limit, 

■(eg) 
NL 



they extract their quantity /a = This quantity is related to our / rcs according to /a = 10/ rcs /9 
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where t k is the time at which k\ ~ hi = k exit the horizon, and we will take t to be some 
late time when the modes with comoving momenta k%, /c 2 and k 3 have exited the horizon. 
In our conventions the two-point function is given by 

^(kx.t^M)) = (27r) 3 5(k 1 + k)|^ o) | 2 . (3.32) 

Recall that 



2 2&l(k) 
= 2n 2 ^^, 3.33 
k 



and remember that A^(fc) is momentum independent in the scale-invariant limit. It acquires 
its momentum dependence from the time dependence of the background. Modes with dif- 
ferent momenta feel a different background as they exit the horizon. In equation (I3.3ip . we 
can thus make the replacement 

1 d -+-/t, (3-34) 



H(tk) dtk dink' 

with the derivative in Ink only acting on the momentum dependence in A^(fc). Restoring 
the k3 dependence inside the 5-function, the consistency relation can then be written in the 
form 



Km Wki, t)K(k 2 , t)K(k 3 , t)) ~ -(2vr) 3 5 3 (k 1 + k 2 + k 3 ) |< } \ 2 \TZ k o) | 2 ^^ • (3-35) 
The logarithmic derivative of the amplitude of scalar fluctuations ( I2.3ip is given by 



dlnA 2 y/W* f(p k \ f2nf\ 1/2 

v ^ J n s — 1 + on s — — sin — with on s = 3o* . (3.36) 



dink f \f J \V2e 

The consistency condition (I3.35j) together with the power spectrum given in Section 2 implies 
that the shape in this limit up to a phase in the trigonometric function should take the form 

g( k ,k,k 3 ) _3V5Fft. ( f ^ 1/2 2^ os ^ ln2 , A ^ . (3 . 37) 



k 2 h 8 \V^J h V / 

This agrees with our result for the resonant shape (I3.17P after setting ki = k 2 = k and taking 
the limit k 3 ^ik. 

Notice that in [13J, the consistency relation f !3.3ip was used to predict the shape of 
resonant non-Gaussianities in the squeezed limit. Here we have derived the resonant shape 
from first principles, and we use the consistency relation as a check of our computation. 



4. Correlations Between Resonant and Other Types of 
Non-Gaussianities 

The comparison of theoretical models of non-Gaussianity with the data is computationally 
very challenging. In light of this difficulty, the authors of [H] have proposed to use a 
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normalized scalar product, or "cosine", to assess to which extent two different 3D primordial 
shapes give rise to a similar 2D signal in the CMB. If two shapes have a large cosine (in 
absolute value), the observational constraints on one can be exported to the other. 

Following this idea, in pT5J, an appropriately defined cosine has been used to classify 
known non-Gaussian models. One class was reserved for models of canonically normalized 
single-field inflation where some feature is present on the top of an otherwise slow-roll flat 
potential. The models considered in this work belong to this class. 

In this section we show that the shape of resonant non-Gaussianity that we have de- 
rived is very different from the shapes of non-Gaussianity that have been constrained by 
data. In particular we compute the correlation (to be defined soon) between resonant non- 
Gaussianity and local, equilateral and orthogonal non-Gaussianity and find that it is always 
less than about 10%. The observational constraints on these models are therefore not useful 
to constrain resonant non-Gaussianity. 



4.1. Scalar product, cosine and shapes 

In this subsection we give the definition of the cosine that we will use in the next subsections 
to compute the correlation between resonant non-Gaussianity and local, equilateral and 
orthogonal shapes. 

Following [T5] , we choose the simplest product that exhibits the same scaling as the 
optimal CMB estimator. The definition is 

F(S, S') = J S(ki, k 2 , h)S'(k h k 2 , k 3 )^ , (4.1) 

where dV = dk\dk 2 dk 3 and 

c / 7 , i \ G(h,k 2 ,k 3 ) 

S{k 1 ,k 2 ,k 3 ) oc — — — . (4.2) 

K\K 2 K 3 

The normalization is irrelevant for our purpose because we will only be interested in the 
normalized scalar product or cosine of two shapes given by 

C(S S') = F(S,S) ,^ 
y/F(S, S)F(S', S') ' 

Because of the rotational and translational symmetries of the background geometry, the 
volume of integration V is three dimensional. The volume three-form dV and the integra- 
tion boundaries are conveniently written in the coordinates {k,a,/3}, which are related to 
{h,k 2 ,k 3 } by [29] 

k=^ = \{h + k 2 + k 3 ) , k 1 = k(l-f3), (4.4) 
k 2 = ^(l + a + /3) , k 3 = ^(l-a + /3) . (4.5) 
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One virtue of this set of coordinates is that for scale-invariant non-Gaussian models, i.e. 
G(k, k, k) oc k 3 or S(k, k, k) oc k°, the integrals over dk cancel between the numerators and 
denominators in f 14 . 3 f) and we are left with a two dimensional integral. However, resonant 
non-Gaussianity is not scale invariant (due to the sine and cosine in e.g. (I4.6P ) and having 
a three-dimensional integration volume is essential to get meaningful results for the cosine 

The integration boundaries in (14. ip require some discussion. In [15], they were chosen to 
be < k < oo, —(1 — j3) < a < (1 — fi) and < (3 < 1. Notice that a and (3 parametrize a 
triangle that is the base of a tetrahedron in the space {ki, k 2 , k^} with the apex at the origin 
and the semi-perimeter k parameterizing the height. With the above integration limits, the 
product ( 14.11) is typically infinite. Depending on the shapes that are being integrated, there 
can be an IR divergence for example where a = ±1 or (3 = 1, i.e. at the vertices of the 
{a, (3} triangle where one of the k t vanishes. Local non-Gaussianity (14 .8p and resonant non- 
Gaussianity (I4.6P show this divergence for squeezed configurations. In addition there are 
generically also UV divergences from the dk integral^!. 

Physically it is clear that for a given experiment, e.g. observations of the CMB or of large 
scale structure (LSS), there is a finite range of momenta {k min , /c max } that can be probed. The 
cosine (I4.3P is useful if it compares two primordial non-Gaussian shapes over the same range 
of momenta that is probed by a chosen class of experiments. Implementing this is slightly 
subtle because of the three-dimensional nature of the primordial non-Gaussian shapes as 
opposed to the two-dimensional nature of the observations (see e.g. [14115] for a discussion). 
The interesting issue of defining a scalar product suitable for non-scale-invariant shapes is 
beyond the scope of this work, so we will limit ourselves to specify some {k m[n , fc max } range 
of integration for ki and check that our results do not qualitatively depend on this choice. 

Let us now turn to the shapes that we will consider. For resonant non-Gaussianity, we 
will work with the linear potential derived from the string theoretic construction. It is clear 
from (I4.3P that the normalization of 5* is irrelevant. We can thus define 

S reB (K h, h) = sin (^J + M cos j g | . (4.6) 

We would like to know the correlation of S res with other shapes that have already been 
compared to and constrained by observations. If the cosine were close to one for some of 
them, we could export their constraints to resonant non-Gaussianity. The best constraints 
form 7-year WMAP data on local, equilateral, and orthogonal non-Gaussianity at 95% CL 
arc [10 

- 10 < / local < 74 , -214 < / cquil < 266 , -410 < f rtho < 6 . (4.7) 



16 As we said, the UV divergence can be neglected in the case of scale-invariant shapes because the dk 
integral simplifies in the cosine (|4.3|) . 
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For optimal limits from 5-year WMAP data see [30|31j . 
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These shapes are defined by 
Siocai(Aa, fa, fa) 

'S'equil(^l) fa, fa) 



«Sortho(fcl) fa, fa) 



£■3 I 1.3 I h.3 
ru-^ 

/C1/C2/C3 

(fei + fa - fa){fa + fa- fa){fa + fa- fa) 



fa fa fa fa 

SS'equilC&l, fa, fa) -2 



fa fa fa 

fa fa fr? fa 



3/cf 



3k\ 



fa fa fa fa fa fa 



Notice that these shapes are factorizable approximations to the results of the theoretical 
calculations (the reader is referred to [31] for further details). They are good approximations 
in the sense that their cosine with the theoretical shapes is very close to one. 



4.2. Numerical results 

We have numerically calculated the cosine ( 14 .3p between the shape of resonant non-Gaussianity 
(I4.6P and local, equilateral and orthogonal shapes in (14. 8p . We have chosen k m - m = 10~ 4 Mpc -1 
and /c max = 10 _1 Mpc -1 for the IR and UV cutoff, respectively. In order to implement the 
IR cutoff in the numerical calculation we have introduced three exponential damping factors 
exp(— k min /ki) with i — 1,2, 3 in the integral (14. ip . The UV cutoff is taken into account 
using the integration region < k < 3k max /2. The results of this brute force approach are 
shown in Figure [31 We have plotted the value of the three cosines as function of the axion 
decay constant /. 

It is clear from the figure that resonant non-Gaussianit y h as a correlation smaller than 
about 10% with the other shapes for any interesting valuer 8 ! of /. The cosines take both 
positive and negative values and get closer to zero as / decreases. Both these features can 
be understood analytically and this is the subject of the next subsection. 



4.3. Semi-analytical results 

To better understand the points in Figure El we have calculated the three cosines semi- 
analytically. We adopted a simplification in implementing the IR cutoff ki > k min : we took 
as integration region — (1 — (3) < a < 1 — (3 and k m i n /k max < f3 < 1 — /c m i n //c max . This cuts off 
not only the vertices of the {«,/?} triangle but also the side of the triangle at (3 = 0. Given 
that none of the shapes we are considering has a divergent contribution along a side (as 
would be e.g. the case for flat shapes obtained in the presence of deviations from a Bunch- 
Davies vacuum [321133] ). the result we get is a good approximation to the one obtained from 
cutting off only the vertices. As can be seen in Figure [31 the results of this subsection nicely 

18 The interesting / range was discussed around (|3.30[) . 
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Figure 3: From the bottom left panel clockwise we plot the cosine between the resonant 
shape (14. 6p . and orthogonal, local and equilateral shapes (|4.8|) . evaluated numerically for 50 
points equally spaced in log/. We also plot the enveloping profiles from the semi-analytical 
calculation to show that numerical and semi-analytical approaches give consistent results. 
In the bottom right panel, numerical results for all shapes are shown together. This makes it 
graphically clear that resonant non-Gaussianity is essentially orthogonal to the other shapes. 



agree with those of the numerical approach described above, in which only the vertices had 
been cut off by the exponential damping factors exp(— fc min /fcj). We have also checked that 
varying the cutoff in /3 away from 10~ 3 by up to an order of magnitude changes the value 
of the cosines by less than a per cent. The cosines plotted in Figure H] as function of / 
are computed using the above integration boundaries for a and /3 and an integration range 

3/Cmin 5: k 5; 3fc max /2 for k. 

The very rapid oscillations are due to the fc-independent phase f and do not have 
a particular physical relevance given that the phase of the oscillations is arbitrary in the 
present model. On the other hand, the enveloping profiles, which are highlighted in Figures [3] 
and HI carry some interesting information. They show that resonant non-Gaussianity has 
a correlation smaller than about 10% with local, equilateral and orthogonal models. The 
observational constraints on the latter are hence of little use in constraining resonant non- 
Gaussianity, the more so the smaller /. The enveloping profiles are well described by a linear 
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Figure 4: From the bottom left panel clockwise we plot the cosine between the resonant 
shape (14. 6p . and orthogonal, local, and equilateral shapes (|4.8|) . as a function of /. In the 
bottom right panel, showing all the three embedding profiles at the same time, makes it 
graphically clear that resonant non-Gaussianity has less than about 10% correlation with 
the other shapes. A linear fit to the enveloping curves is given in (14. 9p . 



fit that gives 

\C(S res ,S local )\ < 1.9/0,, (4.9) 

|C(^e S ,^equil)| < 2.3/0, (4.10) 

\C(S res ,S OItho )\ < 1.3/0,. (4.11) 

This dependence can be understood as follows. Let us consider the cosine C(S res ,S) for 
some slowly varying shape S. There are three scalar products appearing in the definition 
(14. 3 p of C(S res , S). Because of the oscillations of S res , the scalar product F(S res , S) can get 
a contribution from at most half a period, because an integration over one or more whole 
periods gives approximately zero. Notice that for any scalar product the largest contribution 
comes from the region around k ~ /c max because the overall scaling of (14. ip is F(S, S') ~ k^ x . 
The periodicity of S res as function of k around k ~ k max is given by 27r/0,fc max . Hence 
F(S res , S) is obtained effectively from integrating over a range of k that is smaller than half 
a period, i.e. vr/0,/c max . On the other hand, the scalar products F(S res , S res ) and F(S,S) 
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have either squared oscillations or no oscillations at all, respectively. This means that there is 
now no cancellation due to the oscillations and F(S res , S res ) and F(S, S) get integrated over 
the whole range of k, i.e. approximately k max . Taking the ratio as in (14.31) . we see that the 
absolute value of the cosine C(S res , S) can be bounded from above by up to a number 

smaller than but of order one. Let us say it in other words. If there were no oscillations, S res 
would be roughly approximated by some linear combination of local and equilateral shapes. 
Then, always ignoring oscillations, C(S res , S) would generically be smaller than but of order 
one. The presence of oscillations gives the leading effect on the numerator F(S res , S) where 
the range of the dk integral is reduced by a factor 7r/0*. Hence, we find again 



The result of this heuristic argument nicely agrees with the fit of the semi-analytic com- 
putation presented in ( 14.91) . Notice that the argument given above applies to the cosine of 
resonant non-Gaussianity with any slowly varying non-Gaussian shape and not just with 
those considered here. 

5. Discussion 

We have studied the primordial bispectrum of scalar perturbations for models whose potential 
possesses small modulations. We do not deny that our work was largely motivated by a class 
of models derived from string theory that are based on axion monodromy in which such 
periodic modulations on top of an otherwise flat potential are a generic feature [81191134"] . 
However, we argue that these are by no means the only models where such oscillations are 
expected to arise. In large field models of inflation, the inflaton potential is required to be 
flat over a range in field space much larger than Mp\. From the point of view of effective 
field theory, a potential that is flat over such a large range seems unnatural unless there 
is an underlying shift symmetry. This makes axions natural candidates for the inflaton 
especially in the context of large field inflation. If the inflaton is an axion, it seems plausible 
that the potential will receive small periodic contributions from non-perturbative effects. 
These periodic contributions might be due to instantons in a gauge sector the axion couples 
to, or, in the context of string theory, they might arise from Euclidean branes or world- 
sheet instantons. String inspired or not, as soon as we use the shift symmetry of axions to 
explain why the inflaton potential is so flat, we should admit the possibility of small periodic 
modulations in the potential which may lead to observational consequences. So if theoretical 
prejudices have to be employed to isolate a handful of shapes of non-Gaussianity that should 
be looked for in the data, resonant non-Gaussianity deserves to be one of them. We hope 
that, even though it is not factorizable, the analytical expression for the shape of resonant 
non-Gaussianity given in this work will make it possible to obtain observational bounds on 
this type of non-Gaussianity. 

The CMB data is compatible with a small logarithmic running of the power spectrum 
of primordial fluctuations. This has made it natural to consider phenomenological types of 



\C(S res ,S)\ <7T/0,. 
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non-Gaussianity that deviate from scale invariance by at most a small logarithmic running. 
On the other hand, the theoretical considerations expressed above suggest that we should 
keep in mind the possibility of a scale dependence that reproduces scale invariance only after 
an adequate average. 

Resonant non-Gaussianity typically comes with oscillations in the two-point function 
which are constrained by observations |9j. Remarkably, keeping 6* fixed, the amplitude of 
modulations in the two- and three-point function scales in opposite directions when varying 
the frequency. This implies that if oscillations were really imprinted on cosmological per- 
turbations during inflation, they could equally well become observationally accessible in the 
two-point function, the three-point function, or in both. This disentanglement of two- and 
three-point functions is a peculiar feature of resonantly produced perturbations. 

In [IB] , an effective field theory for the fluctuations around a quasi de Sitter background 
was constructed for the case of a single field. As is well known, for a single scalar field 
coupled to gravity, one can fix a gauge in which the fluctuations in the scalar field vanish, or 
are eaten up by the metric. As pointed out in [IS] , this gauge is very much like unitary gauge 
in a spontaneously broken gauge theory where the Goldstone has become the longitudinal 
mode of the gauge field. Also very much like in the gauge theory example, the longitudinal 
mode, or the Goldstone boson dominates the dynamics at high energies so that the non- 
linear theory of the Goldstone contains all the information about the system at sufficiently 
high energies. The theory of the Goldstone, even though non-linear, is easier to study and in 
particular makes relations between different operators transparent that would otherwise be 
obscure. The effective field theory of inflation as presented in [16] or [IT] is the tool of choice 
if one is interested in high energy corrections, i.e. higher derivative corrections. However, 
this high energy limit is equivalent to a limit in which all slow-roll parameters are taken 
to zero. Since the effect of resonant non-Gaussianity arises from a term in the interaction 
Hamiltonian that is higher order in the slow-roll expansion, it should not be surprising that 
it is not captured by the simplest version of the effective field theory of inflation. These 
terms could of course be kept [35J, but essentially at the cost of turning the effective field 
theory of inflation back into the system of a single scalar field coupled to gravity that we 
have studied here, written in a slightly different notation. 

Finally, let us conclude with a couple of interesting directions for future research. We 
have focused our efforts on the three-point function in this work because it is the obvious 
observable to look for when looking for a departure from Gaussianity. The four-point function 
may also be of phenomenological interest in these models, and it can be calculated by the 
same methods presented here. 

Our calculations have shown that the three-point function of primordial curvature per- 
turbations may be large in models with periodically modulated potentials. For a comparison 
with the data, it still remains to calculate the prediction of the model for the two-dimensional 
image of the cosmic microwave background. We have seen that our shape is not factorizable. 
This makes a direct numerical evaluation too time consuming. However, it may be possible 
to make analytic progress in this direction. 
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Constraints on resonant non-Gaussianity could also arise from its effect on large scale 
structures. A very preliminary analysis shows that the effect on the halo bias discussed in 
[36] is relatively modest because it comes from the signal in squeezed configuration which 
is suppressed in our model by the small factor f / yJ2e*. It would be interesting to consider 
other large scale structure observables. 



Acknowledgments 

It is a pleasure to thank Richard Easther, Eiichiro Komatsu, Michele Liguori, Eugene Lim, 
Liam McAllister, Emiliano Sefusatti, and Gang Xu for many useful comments and discus- 
sions. The work of R.F. has been partially supported by the National Science Foundation 
under Grant No. NSF-PHY-0747868 and the Department of Energy under Grant No. DE- 
FG02-92ER-40704. The research of E.P. was supported in part by the National Science 
Foundation through grant NSF-PHY-0757868. 



A. Background Solution for the General Potential 

In this appendix, we derive the background solution and the slow-roll parameters for a 
potential of the form 



y(0) = V o (0)+A 4 cos(JJ , (Al) 

where we assume Vo(0) to be a smooth featureless potential that admits slow-roll inflation. 
For this potential, the equation of motion for the inflaton becomes 

+ 3H( p + v'^) = sin ff \ . (A.2) 

Since oscillations have not yet been observed, we know that the modulations must at least 
be small around the values the inflaton takes when the observable modes exit the horizon. 
This makes it natural to treat the oscillatory part as a perturbation and expand the field as 
= 0o + 0i + ' ' ' • Here 0o is the solution in the absence of modulations, 0i is linear in the 
modulations, and so on. We will assume that 0o is well approximated by the slow-roll result 



K(<f>o) 



(A.3) 



It is convenient to write the equation of motion for <f>i using O as independent variable 
instead of time. To leading order in the potential slow-roll parameters ey and r]v (defined 
in equation ( IB. II) ). this equation can be written a£§ 

/ // 3 /3 3r] Vo \ 3A 4 . /0 O \ 

* - t^, - U - & J * = vi« sm It ) ■ (AA) 



19 This is only true provided \4>i/ f\ <S 1) but we shall see that this is the case. 
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where the slow-roll parameters are functions of 0o and the prime indicates derivatives with 
respect to M Provided we are interested in a motion of 0o that is large compared to the 
decay constant / but small compared to the value of the inflaton when the modes that we 
observe in the CMB exit the horizon, we can replace 0o by everywhere except in the 
argument of the sine. To leading order in slow-roll parameters and assuming / <C v^e*, 
where we use the notation e* = ey (</>*), the result of the equation obtained making the 
above substitution is 



3A 4 / 



2e,V r ( 



sm 



Mt) 
f 



(A.5) 



In analogy to the linear case (11 .11) . it may be convenient to introduce a parameter that 
measures the monotonicity of the potential near the pivot scale. We define it as 



(A.6) 



where the asterisk indicates that, for all potentials except the linear potential, it depends on 
the pivot scale. Using the definition of the potential slow- roll parameters (IB.ll) . the solution 
for the scalar field to linear order in the oscillations can then be written aa 20 l 



<t>{t) = Mt) 



SbJ 2 
\/2eT 



sm 



Mt) 
f 



(A.7) 



with 4>o(t) obtained from integration of equation (1A.3I) . 

With this solution for the background field, it is now straightforward to calculate the 
Hubble slow-roll parameters (12. 2p to leading order in the oscillations. To do this, it is helpful 
to remember the exact relation H = —<f) 2 /2. One finds that to leading order in the slow-roll 
parameters ey and r)y as well as to leading order in the oscillations, the Hubble slow-roll 
parameters (12. 2 p are 



e = e* — 36* /a/^ cos 



36* sin 



Mt) 
f 

Mt) 
f 



where 77* = r]v ((p*) ■ From these, one obtains 



(A.8) 
(A.9) 

/ ~v/;- (A ' 10) 

As long as we are interested only in small neighborhoods of 0*, these results can be 
immediately translated to potentials with multiplicative corrections of the form 



Si 
H 



36*\/2e* / Mt) 
■ cos 1 



V{4>) = V (<f>) 



1 + A cos 



./ 



(All) 



20 Notice that \<f>i/f\< , 3Kf . < 1 as was needed for our derivation of equation (|A.4[) . 

y/2ev (4>*) 
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by replacing A 4 — > AV r o (0*), or equivalently b*fy/2e^ — > A. One finds that the solution for 
the background scalar field is given by 

<P(t) = Mt)- 3 4 L ^(^ 1 ) , (A.12) 



2e* V / 

leading to the slow-roll parameters 



e = e* -3Acos^^ , (A.13) 
3A , {Mt)\ ( a 



* _ 3A /^*) , (A15) 



and finally 

H f 2 ^\ f 

What remains is to find the dependence of </>o on X = —Kt, or equivalently on the 
conformal time which should be obtained upon integration of its equation of motion. In 
the approximation we have been using in the derivation of the solution for the background 
scalar field, i. e. neglecting the time dependence of e and working to leading order in slow-roll 
parameters, the relation is very simple. The equation of motion then takes the form 

d4> ° V2T*. (A.16) 



dln(-r) 

This can immediately be integrated, and with appropriate choice of initial conditions, one 
finds 

M T ) = <P* + v 7 ^ In t/t* . (A.17) 

Here r* is the conformal time at which the pivot scale exits the horizon and 0* and e* 
are the values of the scalar field and of the slow-roll parameter at that time. Multiplying 
numerator and denominator inside the logarithm by k*K, making use of — k*r* = 1, and 
using X = —Kt, this can be written as 

M x ) = <pK + V2e~* In X, with <p K = 0, - v^hi K/K , (A.18) 

where 4>k is the value of the field when the mode with comoving momentum K exits the 
horizon. We conclude that in this approximation 

5l 36 ^cosf^ + 4^1nxV (A.19) 



H f \ f f 
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B. Slow-roll Parameters 



As a courtesy to the reader, in this appendix we summarize various definitions for the choice 
of slow-roll parameters and provide the exact (i.e. valid beyond the slow-roll approximation) 
relations between them. The potential slow-roll parameters are defined by 

1 fV'\ 2 V" 

e " s 2\v) • " >S T- (B1) 

They are very easy to calculate because one does not need to solve for the actual dynamics 
of the system in order to obtain them. The smallness of ty and rjv tells us that there exists 
a regime (which is often an attractor) in which the system evolves slowly. On the other hand 
the potential slow-roll parameters do not carry any informations about the actual dynamics 
(for example the inflaton could be moving fast over a flat region of the potential). They are 
thus of little help in order to assess the validity of the slow-roll approximation for a certain 
dynamics. More useful parameters are those that describe the time evolution, often referred 
to as Hubble slow-roll parameters. While the first slow-roll parameter is always taken to be 

z» = -Jpi ( B -2) 

for the second there are several different conventions: 

_ e H dine 
m = -77 = -7; — , (B.3) 
en a ma 

„. = JL-^ (B4) 

H Id In if ,„ , 

2HH 2 dlna v ' 

where a is the scale factor and a dot denotes a time derivative. Notice that the slow-roll 
parameters are always dimensionless. A third slow-roll parameter is sometimes introduced 

and it is useful in the discussion of perturbations. Exact relations to convert the various 
slow-roll parameters are 

3-^V_. /3 + W2-e 



VV = 3-6 " 3-6 (R8) 
t)h = 25 + 2e. (B.10) 
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